THE MAGNETO-GEOSTROPHIC EQUATIONS: A SURVEY 
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Abstract. We discuss recent results obtained by the autliors in [FVlIa, FVl lb, FV12, FRV12], regarding 
the analysis of the magneto-geostrophic equation: a model proposed by Moffatt and Loper [ML94, Mof08] to 
study the geodynamo and turbulence in the Earth's fluid core. We conclude this review by indicating some open 
problems around the MG equation, that remain to be addressed in the future. 



Contents 

1 . Introduction 1 

2. Diffusive MG equation 5 

2.1. Global regularity of singular drift-diffusion equations 5 

2.2. Lyapunov stability of steady states 7 

3. Non-diffusive MG equation 11 

3.1. On the (lack of) well-posedness for the MG equation in Sobolev spaces 1 1 

3.2. Hadamard ill-posedness for the non-diffusive MG equation 12 

4. Fractionally diffusive MG equation and well-prepared initial data 14 

4.1. Regime of well-posedness in Sobolev spaces 15 

4.2. Regime of ill-posedness in Sobolev spaces 15 

4.3. The threshold case 7 = 1/2 16 

4.4. Well-prepared initial data 17 

5. Conclusions 18 
Acknowledgements 19 
References 19 



1. Introduction 

Physicists have long realized the importance of the Earth's magnetic field and that this field originates in 
the Earth's fluid core. The geodynamo is the process by which the rotating, convecting, electrically con- 
ducting molten iron in the Earth's fluid core maintains the geomagnetic field against ohmic decay. The 
convective processes in the core that produce the velocity fields required for dynamo action are a combina- 
tion of thermal and compositional convection. The full dynamo problem requires the examination of the full 
three dimensional partial differential equations governing convective, incompressible magnetohydrodyam- 
ics (MHD). In the past decades computer models have been used to simulate the actual geodynamo. See, for 
example, the review articles by Roberts et al on modeling the geodynamo [BR95, GR97, RGOO, RGOl] and 
references therein. However, as Glatzmaier et al [GOC04] remark, no three dimensional dynamo model has 
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been run at the spatial resolution required to simulate the broad spectrum of turbulence which surely exists 
in the Earth's fluid core. 

Cunent computers and numerical methods require the imposition of diffusivities that are several orders of 
magnitude larger than those which are realistic. It is therefore reasonable to attempt to gain some insight into 
the geodynamo by considering a reduction of the full MHD equations to a system that is more tractable, but 
one that retains many of the essential features of the problem. The magnetogeostrophic equation proposed 
by Moffatt and Loper [ML94, Mof08] is one such model. The physical postulates of this model are the 
following: slow cooling of the Earth leads to slow solidification of the liquid metal core onto the solid 
inner core and releases latent heat of solidification that drives compositional convection in the fluid core. 
The arguments for the appropriate ranges of the characteristic length, velocity, and perturbation density are 
based on these physical postulates. 

We first present the full coupled three-dimensional MHD equations for the evolution of the velocity vector 
u{x,t), the magnetic field vector B{x,t) and the buoyancy field 9{x,t) in the Boussinesq approximation 
and written in the frame of reference rotating with angular velocity 17. The physical forces governing this 
system are Coriolis forces, Lorentz forces, and gravity. Following the notation of Moffatt and Loper [ML94] 
we write the equations in terms of dimensionless variables. The orders of magnitude of the resulting non- 
dimensional parameters are motivated by the physical postulates of the Moffatt and Loper model: 

[Ro (dtu + u • Vn) + 63 X n] = -VP + ea • V6 + R„,b ■ Vb + N^Oe^ + eAu (1.1) 

Rm[dtb + u-Vb-b-\/u\= e2-Vu + Ah (1.2) 

dte + u-V9 = e^A6 + S (1.3) 

V-u = 0,V-6 = 0. (1.4) 

Here S{x, t) is a given smooth function that represents the forcing of the system. 

The mathematical statement of the geodynamo problem asks whether there are initial data for the MHD 
system for which the evolution of the perturbation of the magnetic field b{x, t) grows for sufficiently long 
time. This can be interpreted as a question of the existence of instabilities of (1.1)-(1.4). 

The notation in (1.1)-(1.4) is the following. The Cartesian unit vectors are denoted by ei, 62, and 63. 
For simplicity, we have assumed that the axis of rotation and the gravity g are aligned in the direction of 
63. We have assumed that the magnetic field B{x, t) consists of the sum of an underlying constant field 
BQe2 and a perturbation b{x, t). Our choice of -60^2 as the underlying magnetic field is consistent with the 
models where the magnetic field is believed to be predominantly toroidal (cf. [ML94]). The dimensionless 
parameters are 

where the typical velocity scale V is chosen to be 

V^lf. (1.6, 

Here ^0 is a typical amplitude of the buoyancy 6 and L is a characteristic length scale. P is the sum of 
the fluid and magnetic pressures, v is the kinematic viscosity, rj is the magnetic diffusivity, and k, is the 
molecular diffusivity of the compositional variation that creates an ambient density p. The usual estimate 
for 7] in the fluid core is Sm^ cm^^. The values of v and n are speculative, but likely to be extremely small. 
For a more detailed discussion of plausible ranges of the physical parameters that are appropriate for the 
geodynamo, we refer the reader to the book of Ghil and Childress [GC87]. 

For the regions in the Earth's fluid core modeled by Moffatt and Loper, it is argued that the dimensionless 
parameters are the following orders of magnitude. The inverse Elsasser number A^^ will be of order unity, 
or somewhat less in the condition of the fluid core, while the Rossby number Ro is of order 10^'^. The value 
of £ is somewhat uncertain since the viscosity 1/ is a matter of guess work, but plausibly e < 10^^ in the 
core conditions. The inverse Peclet number is also likely to be extremely small, perhaps of the order of 
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10~^. According to Moffatt and Loper [ML94], the magnetic Reynolds number Rm is relatively small in 
the range of L considered. With these estimates of the sizes of the parameters the dominant balance of the 
leading order terms in equations (1.1) and (1.2) give the following reduced system: 

Ar2e3 X M = -VP + 62 • V6 + iV^6'e3 (1.7) 

= 62 • Vu + Aft (1.8) 

along with the incompressibility condition 

V-u = 0,V-6 = 0. (1.9) 

The linear system of equations (1.7)-(1.9) determine the differential operators that relate the vector fields 
u{x, t) and h{x, t) with the scalar buoyancy 6{x, t). These operators encode the vestiges of the physics in 
the problem, namely Coriolis force, Lorentz force, and gravity. Vector manipulations of (1.7)-(1.9) give the 
expression 

[N^{ez ■ V)2A + (62 • V)^] u = iV^(63 • V)A(63 X Ve) + N^{e2 ■ VfV x (63 x V^). (1.10) 

The sole remaining nonlinearity in the system comes from the coupling of (1.10) with the equation (1.4) for 
the time evolution of the 6{x,t). We call this nonlinear active scalar equation the magnetogeostrophic (MG) 
equation. It is the subject of this article. 

We study the properties and behavior of the active scalar equation in T'^ 

dte + u-Ve = e^M + S (1.11) 

u = M[6] (1.12) 

via an examination of the Fourier multiplier operator M obtained from (1. 10), which relates u and 9. More 
precisely, 

Uj = Mj[e] = {MjOy 

for j G {1, 2, 3}, and the explicit expressions for the components of M as functions of the Fourier variable 

k = (fei, A;2, A;3) G '^l ■= \ {ks = 0} are 

Mi{k) = {N'^k2k3\k\'^ - N^kik^) (1.13) 

M^{k) = {-N^kik:i\k\^ - N^k^k^) D'^ (1.14) 

M^{k) = {N^klkl + N^ki) D-^ (1.15) 

where 

D = N^\k\^kl + ki. (1.16) 

On {A;3 = 0} we define Mj{k) = 0, for all j € {1, 2, 3}, since for self-consistency of the model we assume 
that 9 and u have zero vertical mean, i.e. J 9dx3 = J udx^ = 0. Several important properties of the Mfs 
are immediately obvious. Firstly, the functions are strongly anisotropic with respect to the dependence on 
the integers ki, k2, and ^3. This is a consequence of the interplay of the three above mentioned physical 
forces. Secondly, the function M3 is an even function of fci, k2, and ^3. Thirdly, we have Mj{k)kj = for 
any kGZ^ (we use the standard summation convention on repeated indices), which is built into the system 
since the velocity field is divergence free 

V • n = 0. 

Lastly, although the symbols Mj are bounded (by a constant that depends only on N) in regions of Fourier 
space where k2 < max{A;2, ^3}, this is not the case in the "curved" regions where k^ = 0{l), and k2 = 
0{\ki\''') with < r < 1/2. In such regions the symbols are unbounded as |fci| — )• 00 with 

|M(A:)| <C|A;|, (1.17) 

for some suitable chosen constant C that is independent of N, and this bound is sharp. Hence the MG 
equation falls in the class of abstract active scalar equations (see (2.1)-(2.2) below) where the relation 
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between the divergence-free velocity field u and the scalar 6 is given as a singular operator of order 1. These 
properties of M make the MG equation interesting and challenging mathematically, as well as having a clear 
physical basis in its derivation from the MHD equations. 




Figure 1. From left to right: plotof Mi(A;i, k2, 1), M2{ki, k2, 1), and Ms{ki, k2, 1), with 
N = 1. The range of the frequencies is ki £ [-10^, 10^] and k2 € [-10^, 10^]. 

In Section 2 we state results proved in [FVlla, FVllb, FV12] for the diffusive MG equation, i.e., the 
case when > 0. We recall that in the fluid core is positive but probably extremely small. In this case 
the MG equation is global well posed. To prove this fact we use the control on H^oo , a quantity that is 
non-increasing along solutions, and which gives a bound on the drift velocity u in the space L'^BMO^^. 
The latter information on the drift is "critical" for the natural scaling of the linear drift-diffusion equation 

dte + vve = M. (1.18) 

Indeed, for v £ LfBMO^^ n L^L^ divergence-free, we proved in [FVlla] that the solution d of (1.18) 
becomes Holder continuous in finite time, by appealing to the classical ideas of De Giorgi [DG57], in the 
spirit of [CVIO]. This result was also established independently in [SSSZ12]. See also [Nas58, Osa87, 
Zha04], the more recent [NU12, SV12, Fill2, IKR12], and references therein, for regularity issues of linear 
parabolic equations with singular coefficients and divergence-free drift. The Holder continuity of solutions 
to the linear equation (1.18) may then be bootstrapped to show that the solutions of the dissipative MG 
equations become smooth in arbitrarily short finite time [FVl lb]. 

Having established the global well-posedness of (1.11)-(1.12) for any > 0, we can proceed to examine 
the stability of equilibria. We show that the existence of an unstable eigenvalue for the linearized equations 
implies nonlinear (Lyapunov) instability of the nonlinear MG equation (see also [FSV97, BGS02, VF03, 
FPS06, FPV09]). We then exhibit the existence of an unstable eigenvalue in a particular example, using the 
method of continued fractions cf. [MS61, FSV97]. The magnitude of this eigenvalue is 0{e'^^) which, in 
the context of the fluid core, is extremely large. Thus the MG equation can be used to model very rapid 
exponential growth in time of 9, and hence, through (1.1 1)-(1.12), of the perturbation velocity and magnetic 
fields, i.e. a dynamo type instability. We emphasize that this "MG-dynamo" instability requires the presence 
of all three forces - Coriolis, Lorentz, and gravity - because the dimensionless parameter N"^ (defined above 
in (1.5)) must be strictly in the range (0, oo). 

In Section 3 we examine the non diffusive MG equation and set 6^ = (cf. [FVllb]). The situation is 
then dramatically different from the case > 0. Without the Laplacian to control the unbounded operator 
M, the MG equation is ill-posed in the sense of Hadamard in Sobolev spaces. Such a result utilizes in an 
essential way both the unboundedness of and the fact that it is an even function of k. The proof again 
uses the existence of an unstable eigenvalue of the linearized equation, however when = it is possible 
to construct such an eigenvalue that is arbitrarily large and it then follows that there is no Lipschitz solution 
map at time t = (see also [Tao06, Ren09, GVDIO, GVNIO, GNIO]). 

In Section 4 we describe results for the fractionally diffusive MG equation (cf. [FRV12]), i.e. replace the 
Laplacian by —(—A)'*', for 7 G (0, 1). This situation is non-physical but mathematically interesting in the 
context of the tools required to treat the fractional Laplacian. We now have a dichotomy across the value 
7 = 1/2. In the range 1/2 < 7 < 1 the equation is locally well-posed, while it is Hadamard Lipschitz 
ill-posed for < 7 < 1/2. At the critical value 7 = 1/2 the problem is globally well-posed for suitably 
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small initial data (and source term), but is ill-posed for sufficiently large initial data. In terms of the non- 
dimensional parameters introduced above in (1.5), this dichotomy can be seen as global well-posedness for 
^ 1, and Lipschitz ill-posedness for <C 1. A further feature of interest is that the anisotropy of the 
symbol can be explored to obtain an improvement in the regularity of the solutions when the initial data 
and source are supported on a plane in Fourier space. For such well-prepared data the local existence and 
uniqueness of solutions can be obtained for all values 7 G (0, 1), and the global existence holds for all initial 
data when 7 E (1/2, 1). 

We note that some of the theorems discussed in this article are valid for certain classes of partial differ- 
ential equations cf. (2.1)-(2.2) below, and in any spatial dimension d > 2. For such theorems, the three 
dimensional MG equation is a particular example of an active scalar equation for which these results hold. 

We conclude the manuscript with Section 5, in which we present a few remaining open problems regard- 
ing the mathematical analysis of the MG equations. 



In this section we discuss the global regularity of solutions to (l.ll)-(l. 12) and the Lyapunov stability of 
steady states, in the case that > 0. 

2.1. Global regularity of singular drift-diffusion equations. In order to show that there is no finite time 
loss of regularity of solutions to (1.1 1)-(1. 12) with > let us first perform a rough derivative count. The 
dissipation is smoothing by two derivatives, while the nonlinearity is also losing two derivatives: the first 
coming from V9, while the second one derivative loss is coming from the constitutive law u = M[9], in 
view of (1. 17). In this sense we may call the problem "critical" (we make this more precise below). 

Since the operator M relating u and 9 is an order one pseudo-differenential operator, one may view the 
MG equation as an example of a more general class of of active scalar equations: 



with Tij a matrix of zero-order pseudo-differential operators, and we use the summation convention on 
repeated indices. We have dropped the smooth source term S from (1.11) since it is not an obstacle to 
the well-posedness theory, have re-scaled the time variable so that the diffusivity coefficient is 1 instead 
of (this is not an issue since we will construct global in time solutions), and the Tij corresponding to 
(1.12) are given by Tij = —di{—A)~^Mj. We call this class of active scalar equation "singular" since 
the drift velocity is given in terms of the advected scalar by a constitutive law which is losing derivatives. 
We refer the interested reader in to the work of Chae, Constantin, Cordoba, Gancedo, and Wu [CCC+12] 
and also [FGSV12, MX 12] for further discussions about the well-posedness theory of singular active scalar 
equations. 

To make the idea of criticality more precise, we note that the active scalar equation (2. l)-(2.2) has a 
natural scaling invariance built into it: if 9{t, x) is a solution, then 9x{t, x) = 9{\^t, Ax) is also a solution, 
with corresponding drift given by u\{t, x) = \u{\^t, Ax) = diTij9x{t, x). Hence L°°(M'^) is the critical 
Lebesgue space with respect to the natural scaling of the equation. Additionally (2.1) is a scalar equation, 
so that LP norms are non-increasing in time for 1 < p < 00, and hence the norm is the "strongest" a 
priori controlled norm along solutions. In this sense we say that (2.1)-(2.2) is a "critical" equation. 

Another active scalar equation derived from the Navier-Stokes equations for a rapidly rotating fluid in a 
thin shell is the surface quasi-geostrophic equation, introduced by Constantin, Majda, and Tabak [CMT94] 
(see also [CCWOl, CC04, Wu05, KislO] and references therein), which in the critically dissipative case 
reads 



2. Diffusive MG equation 



dt9 + (u • V)9 = A9 

Uj = diTij9, djUj = 



(2.1) 
(2.2) 



dt9 + (n • V)9 
u = V^(-A)- 



= -(-A)^/20, 

^/^9 = {R29,-Ri9) 



(2.3) 



(2.4) 
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and Ri represents the i^^ Riesz transform. There are strong analogies between the systems (2.3)-(2.4) and 
(2. l)-(2.2). The velocity field is divergence free in both cases, and the L°° norm is critical with respect 
to the natural scaling of the equation. The main difference is that in for SQG the dissipative operator is 
non-local, but the drift velocity is not more singular than the advected scalar. 

Inspired by the De Giorgi-type proof of global regularity for the critically dissipative SQG equation of 
Cafarelh and Vasseur [CVIO] (see also [KNV07, KN09, CV12] for different proofs), to show that solutions 
of (2.1)-(2.2) do not blow up in finite time, we first look at the associated linear problem, where the drift 
velocity is given and obeys the bounds that are a priori available for weak solutions of the nonlinear problem: 
u E L°°(0,oo;SAfO-i) nL2(0,oo;L2). In [FVUa] we use a De Giorgi iteration technique to prove that 
a weak solutions of the linear drift-diffusion equation 

dtO + {v ■ V)9 = AO (2.5) 

V-t; = 0, Vj = diVij, Vij e L°^{0,oo;BMO)nL^{0,oo;H^) (2.6) 

instantly become Holder continuous. More precisely, the result is: 

Theorem 2.1 ([FVUa]). Let the initial data Oq £ and 6 € L°°{0, oo; L^) n L^(0, oo; H'^) be a weak 
solution of (2.5), where the drift velocity v is given by {2.6) for a suitable matrix Vij. Then, for any to, there 
exists a > such that 9 G C"/^(to, oo; C"). 

The above result was obtained independently by Sergin, Silvestre, Sverak, and Zlatos in [SSSZ12]. In 
view of the constitutive (2.2) a direct corollary of Therem 2.1 is that: 

Corollary 2.2. Weak solutions of (2.1)-(2.2) are Holder continuous in positive time, i.e. they obey 6 G 
L°° (to, oo; C°')for any to > 0, and a suitable a > 0. 

In view of the above discussion about scaling, Holder regularity is sub-critical information for (2.1)-(2.2) 
and one may bootstrap to prove that the solution is classical: 

Theorem 2.3 ([FV12]). Let G L°°{tQ, oo; C") be a weak solution of the initial value problem associated 
to (2.1)-(2.2), with a G (0, 1). Then 6 G oo; C^+^),forany ti > to, for some 6 G (0, 1). 

The proof of Theorem 2.3 relies on a suitable bootstrap procedure in Besov spaces, in the spirit of [CW08] . 
Repeating the argument for the derivative of 6 one may further obtain that the solution of (2.1)-(2.2) is in 
fact C°° smooth in finite time. 

We conclude this subsection with a very brief discussion of the proof of Theorem 2.1. We refer the reader 
to [FVUa] for all details. Along the lines of Caffarelli and Vasseur [CVIO, Theorem 3] for the critical 
SQG equation, the first step consists of showing that a weak solution is bounded in the spacial variable 
(see [FVUa, Appendix A] for a construction of weak solutions which obey the energy inequality). This 
is achieved by proving that the function {0 — = max{^ — h, 0} has zero energy if h is chosen large 
enough. The key ingredient is the inequality 

dt{e - h)+ - A{9 - h)+ + {v ■ V){e -h)+<0 

which can be shown to hold for weak solutions. Using that v is divergence free one obtains the energy 
inequality 

/ \{e{t2r) -h)+\'^dx + 2 r ( W^\V{e -h)+\'^dxdt< I \{e{ti,-)-h)+\'^dx, (2.7) 
Jm.d Jti J JR<* 

for all /i > and < ti < t2 < oo. For to > 0, and > to be chosen sufficiently large, letting 
tn = to- to/2'', hn = H- i7/2", and 

/" roc r 

Cn = sup \{e{t,-) - hn)+\'^dx + 2 / \V {6 - hn)\'^ dxdt, 

t>t„ jRd Jt J^d 
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it follows from (2.7) that 



C 



n{l+4/d) 1+2/d 



for some constant C > 0. Letting H = H{\\9o\\j^2^^d^,tQ^^) be sufficiently large we thus obtain implies that 
— )• exponentially as n — )■ cxd, and hence 9{tQ,-) < H. 

Having established the boundedness of the weak solutions, the next step is to prove Holder regularity, 
which requires a more delicate argument. The main obstruction to applying the classical parabolic De 
Giorgi estimates via an -based Caccioppoli inequality (1 < p < oo), is that Vj = diVij is the derivative of 
a BMO function. By essentially using the divergence-free nature of v, and appealing to the John-Nirenberg 
inequality, we may however prove a suitable local energy inequality: for any < r < i? and /i G M, we 
have 

ll(^-/^)+llW(Q.) + l|V(^-/.)+lli.^(Q.) < (2-8) 

where C = C{d, \\Vij\\LfBMOx) ^ fixed positive constant. Here we have denoted by Qp the parabolic 
cylinder [to — P^, ^o] x Bp{xo) for /? > and (to, xq) G (0, oo) x M''. Moreover, estimate (2.8) also holds 
with 9 replaced by —9. The bound (2.8) may be in turn used, using a suitable iteration argument to obtain 
an improvement of oscillation result. Assume that ho < supg^^^ 9, where ro > is arbitrary. Then we have 

sup.<.„ + c(H^>"°'"'^"l'"^""V''Lp.-J (2.9, 



Qro/2 V / \<3 



for some positive constant C = C{d, \\Vij\\L^BMOx)- We note that the prefactor of supg^^^ 9 — ho of the 
right side of (2.9) does not scale with r^. Besides estimate (2.9), the main ingredient in the proof of Holder 
continuity is the following estimate on the growth of level sets of the solution. There exist kq G (0, 1), 
^0 = <^o('«Oi \\yij\\LfBMOx)^ ^iid no > 2 such that for ti, i? > if we know that 

> h}r\Br\ < l-\Br\, 



then for all t2 G [ti, ti + 5or ] we have 



|{0(t2,-) > ^}nB/j| < (2.10) 

where r = k^R, M = SMV(^t^^t^+s^K^)^Ba^^ = '^'^^{t^,t^+5oR^)xBn&^ h = {M + m)/2, and H = 
M — {M — m)l 2"" . To prove the above estimate one uses a suitable version of the local parabolic energy 
inequality. 

Bounds (2.9) and (2.10) are the main ingredients in the De Giorgi iteration procedure, which at the end 
of the day yields the improvement in oscillation 

I sup0-infe 1 < 7 I sup0-inf^ 1 (2.11) 

for some 7 G (0, 1), where Qi = [ti, ti + ^o''^] x B4, Q2 = [ti, ti + SoR"^] x B^, ti, R > are arbitrary, 
r = kqR, and 60, kq G (0, 1) is as above. Since 7 < 1, is then classical to show that (2.11) implies the 
solution 9 is Holder continuous, with exponent that may be computed explicitly. 

2.2. Lyapunov stability of steady states. In this subsection we consider the stability of steady states to 
(1.11)-(1.12). First we prove for steady states of the the more general class of equations (1.11), with drift 
velocity given by (2.2), that linear implies nonlinear instability. We then turn to the specific MG equations 
(1.11)-(1.12), for which we construct an explicit unstable eigenvalue of the linear operator. 
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2.2.1. Linear implies nonlinear instability. Let ©o G C°° be a steady solution of (l.ll)-(2.2), i.e. 

f/o • VBo = e.AGo + 5, Uoj = diTijQ^, V • C/q = 0, (2.12) 

where S E C°° is a time independent force, > 0, and Tij are zero order pseudo differential operators. 
Writing the solution of (l.ll)-(2.2) as + 0o we obtain that the perturbation temperature 9 obeys 

dt9 = Le + N0, (2.13) 

where the dissipative linear operator L is defined as 

Le = -Uo-Ve -u-Veo + e^A6 (2.14) 

with Uj = diTijO, and the nonlinear operator N is given by 

Ne = -u-ve = -v -{ue). (2.15) 

Recall the following notion of Lyapunov stability: 

Definition 2.4 ( [FS V97]). Let {X, Z) be a pair of Banach spaces. A steady state ©o is called {X, Z) 
nonlinearly stable if for any p > 0, there exists p > so that if9{-,0) € X and \\0{-,0)\\z < p, then 

(i) there exists a global in time solution to the initialvalue problem (2.13)-(2.15) with G C([0, oo); X); 

(ii) and we have the bound \\9 {-,1)11 z < pfora.e. t G [0,cx)). 

An equilibrium ©o that is not stable in the above sense is called Lyapunov unstable. 

In the above definition, the Banach space X is the space where a local existence theorem for the nonlinear 
equations is available (e.g. X = H'^), while Z is the space where the spectrum of the linear operator is 
analyzed, and where the instability is measured (e.g. Z = L^). The following result states that the linear 
Lyapunov instability of a steady state ©q implies the nonlinear one. 

Theorem 2.5 ([FVl lb]). Let ©o be a smooth, mean zero solution of (2. 12). If the associated linear operator 
L, as defined in (2.14), has spectrum in the unstable region, then the steady state is {H^,L'^) Lyapunov 
nonlinearly unstable, for arbitrary s > d/2. 

The proof relies on a bootstrap argument that is by now well-known for fluids equations [FS V97, BGS02, 
VF03, FPS06, FPV09], and an a priori time independent bound on the norm of the solution, which we 
obtain from the De Giorgi C" estimate. 

We only briefly sketch the argument given in full detail in [FVl lb], to which refer the interested reader. 
First, we fix a linearly unstable smooth eigenfunction (p, with eigenvalue of maximal real part A. In order to 
prove Theorem 2.5, for e G (0, we consider solutions of 

dt9' = Le' + N{9'), e'\t=o = £'P (2.16) 

where L, N, and </> are as above, and show that there exits a constant > such that for any e G (0, e*] 
we may find a time > with 

W{T,)\\l2>C,. (2.17) 

By Definition 2.4, the bound (2.17) shows that the trivial solution 6* = of (2.13)-(2.15) is 
Lyapunov unstable: indeed, an 0{e) perturbation leads in finite time to an 0(1) size solution, in finite time. 

In order to prove (2.17) we use the following classical idea: the initial data being an eigenfunction of L, 
will grow exponentially in time, while the nonlinear term may be suitably controlled iso that it cannot cancel 
this exponential growth. To make this idea concrete, it is convenient to define Ls = L — {X + 6)1, where 
< (5 < A/8. Due to the shift, the resolvent of Ls contains the full right-half of the complex plane, and 
being a bounded perturbation of the Laplacian we also have that Ls generates an analytic semigroup over 
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L^, and ^ is smoothing by two derivatives. Then from the Duhamel formula and the interpolation bound 

11^5^/^(0^)11^2 < C7||r||^/2^||Ar||^//, we obtain the estimate 



< |ke^VlU2+ / e(^+^)(*-^)||L,'/V*(*-^)||£(i2)||L-'/V(r(s))||,.2(i. 

JO 

< C^ee^' + C fe(^+')(^-^) ^ ^^j e^s)\\%'\\Ae^s)f/,'ds. (2.18) 

Jo [t-s)'/^ 

The key step which allows us to close the estimates is proving the existence of a constant C* = C* (©o, S) 
such that 

||Ar||i2 < C* < oo (2.19) 

for all sufficiently small e. The proof of estimate (2. 19) relies on the C" bound which may be obtained from 
Theorem 2. 1 and a suitable bootstrap procedure to reach the regularity level (see [FV lib, Appendix A]). 
Lastly, let i? > := ||0||l2 and T = T{R, e) be the maximal time such that 

L2 < ei?e^* for all < t < T. 



0^{t)\\L^ < C^ee^* + Ci (eRe^^Y^^ (2.20) 



From (2.18) and (2.19) one may then obtain 

^9/8 

and then in turn that 

T>Te = \~^\og{C2/e). 

A lower bound on ||0'^(f)||j^2, similar to (2.20) but with the linear term minus the 9/8 power term, may be 
obtained from the Duhamel formula. In turn, this implies that ||(?^(r£)||^2 > C2{2C^ — i?) =: C* > and 
concludes the proof of Theorem 2.5. 

2.2.2. The existence of unstable eigenvalues for the linearized MG operator. In order to apply Theorem 2.5 
to the (nonlinear) MG equations ( 1 . 1 1)-( 1 . 12), we need to find a steady state ©o for which the corresponding 
linear operator L defined in (2. 14), has unstable spectrum. 

To achieve this we use an idea that dates back at least to [MS61] for Navier-Stokes and [FSV97] for 
Euler. First, note that ©o = -^(2:3), with F smooth, is automatically a steady state, with associated velocity 
Uq = 0, if we assume that the source term is given by 5 = —ei^F"{x^). We thus consider the source term 
to equal S = St^m? sin(mx3), for which we have a steady state solution of (1.1 1) given by 

©0 = sin(mx3), Uq = 0. (2.21) 

For the above choice of the steady states, the equation dtO = L6 becomes 

dt6{x,t) = L6{x,t) = —mu-i{x.,t)cos{mxz) + et^9{x,t), u^ = M^[0], (2.22) 

where M3 is the Fourier multiplier defined earlier in (1.15)-(1.16). Next, we make the ansatz that the 
temperature 9 which solves (2.22) has the form 

d{x, t) = e*^* sin(A;iXi) sin(A;2X2) Cp sin(mpx3) (2.23) 

p>i 

for some fixed ki, k2 in to be chosen precisely later. Using (1.15)-(1.16) we may compute from the 
third component of the velocity field u associated with 9 given by (2.23), insert this into (2.22) and obtain 
that a and {cp}p>i must solve 

fjp Cp smimpxs) + — (sin(m(p + 1)2;3) + sm(m(p — l)x3)) = (2.24) 
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where we have denoted 

o-p = a + {kl + kl + {mpf) (2.25) 

and 

— = mN^ kUkf + kj) 

ap N^{mp)^{kf + k^ + {mpf) + kl ^ ' 

with N being the non-dimensional constant in (1.5). An essential feature of the coefficients ap is that they 
form an increasing sequence which grows rapidly, as when p — )■ oo. Equation (2.24) gives a recurrence 
relation for the sequence {cp}p>i (see e.g. [FVl lb, FRV12]): 

(TpCp + + = 0, for p > 2, (2.27) 

(Jici + — = 0, for p = 1. (2.28) 

We emphasize that given any a > (which then uniquely defines ap for all p > 1) and given any ci > 0, 
one can solve for all {cp}p>2 using the recursion relations (2.27)-(2.28). The key is however that only 
for suitable values of a do the Cp vanish sufficiently fast so that 9 defined by (2.23) is sufficiently smooth. 
Solving for such a amounts to solving the infinite continued fraction equation 

aiai = ^—^ . (2.29) 

0-404-... 

Solving (2.29) is the main difficulty in this argument, since we need not just show there is a solution, but also 
give bounds on it in terms of the Op. We omit these details and summarize the results obtained in [FVl lb, 
Lemma 2.6], [FRV12, Lemma 4.2], and [FGSV12, Theorem 2.2] as follows. 

Lemma 2.6 (Bound of solution to continued fraction equation). Let ap, ap be defined by (2.25) and 
(2.26) respectively, and let ci = ai. There exits a real eigenvalue o"* = (t*(A;i, k2, m, N, e^) > which 
solves (2.29), and we have the estimate 

1 2 

< (^1(^2 < • (2.30) 

aia2 a\a2 

Moreover, the corresponding sequence {cp}p>i decays exponentially in p. 

A direct consequence of estimate (2.30) is that the eigenvalue cj* obeys the upper and lower bounds 

For N, 6^,771 fixed, we hence obtain from (2.32) that there exists at most finitely many frequency pairs 
(fci, /C2) € 1^1 for which cr* > 0. Moreover, one may show that in the parameter range of the geodynamo, 
i.e. < Ek <^ 1, for suitable {ki, ^2) the left side of (2.31) is positive (see (2.33) below), thereby proving: 

Theorem 2.7 ([FVl lb]). Associated to the steady state Bq defined in (2.21), the linearized MG operator L 
defined in (2.22), has unstable spectrum. 

Lastly, note that one may optimize in {ki,k2), and compute a lower bound for the largest eigenvalue of 
L in terms of the physical parameters N, £«, and of m. More precisely, we have that 

9 AAT4 2rJ\ ,2^y 2^ ^ I.A " + ^1 + ^ ^ " ^^.33) 
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with the lower bound being attained e.g. for (of course, we need to take integer part) 

1 NJm 

and k2 = , for < <C 1. 



In particular, we note that as — ^ the lower bound given by the right hand side of (2.33) blows up (this 
will be important in Sections 3 and 4). Moreover, for the above choice of ki and k2, we have kik^'^ — )■ 
l/{2N'^m) as e^. — ?■ 0. This is in precise accordance to the fact that the absolute value of the Fourier 
multiplier associated to M grows at the fastest rate when \k2\ = O^kil^^"^) as \ki\ — oo. See (1.17). 

Combining the results in this subsection with Theorem 2.5, we have thus constructed an example of 
solutions to the MG equations which displays strong (exponential) growth in 6 and hence in the velocity 
field u by (1.12) and also in the magnetic field b{x, t) through (1.8). The rate of this growth can be shown 
to be proportional to exjp{t/2^eK,) which is consistent with the dynamo instability scenario. 

3. Non-diffusive MG equation 

As discussed in the introduction, in the fluid core the value of e^, given by the quotient of k the molecular 
diffusivity coefficient, and L V, the product of the typical length scale and typical velocity in the fluid, is 
very small. 

Therefore, from the physical point of view it is natural to study the non-diffusive (unforced) MG equation 

dtO + {u-V)e = (3.1) 
u = M[0], V-n = (3.2) 

or the coiTcsponding non-diffusive singular active scalar equation when (3.2) is replaced by (2.2), i.e. uj = 
diTijO. As we shall see below, the system (3.1)-(3.2) is ill-posed in Sobolev spaces, in the sense that there 
is no operator which is Lipschitz in time at t = 0. 

3.1. On the (lack of) well-posedness for the MG equation in Sobolev spaces. Before we discuss the 
issue of ill-posedness, let us briefly discuss why the evenness of the operator M breaks the classical proof 
of local existence in Sobolev spaces for the class of singular active scalar equation (3. l)-(3.2). The standard 
energy-approach to local existence in would entail applying the (— A)*/^ to (3. 1), and then take inner 
products with {—Ay/'^O. Morally speaking (this is rigorous when s/2 is an integer in view of the Leibniz 
rule) we have only one term which seem to prevent closing the estimates at the level: 

T = j (-A)^/\ • ve {-Ay/'^e dx = j Mj{-Ay/^e djO {-Ay/^9 dx (3.3) 

Indeed, the other term where more than s derivatives lie on a sing le function is / u-V{-Ay/'^e{-Ay/'^edx, 
but this term vanishes upon integrating by parts since V • li = 0. The only hope to treat the term T would 
be to discover a commutator structure. However, since M is not anti-symmetric, i.e. even in Fourier space, 
we cannot write T = —S, where 

s = j Mj (^{-Ay^^e djOj {-Ay/^e dx = T+ j [Mj,dje]{-Ay^^e {-Ay/^edx = t + c (3.4) 

and [•, •] is the usual commutator bracket. Instead for the MG, we have T = S, and we cannot explore the 
fact that C may be bounded in terms of (if s > 2). This is the main reason why we do not seem to be 

able to close estimates at the Sobolev level. 

On the other hand, it was observed in [CCC^ 12], see also [FVl lb], that if we replace u = M[9] (cf. (3.2)) 
by a certain class of singular constitutive laws Uj = diTijO (cf. (2.2)), and the matrix Tij is such that 

{Tijf,g) = {f,Tijg) 

for any f,g and i,j (here (•, •) is the usual inner product), then the corresponding singular active scalar 
equation (3.1)-(2.2) is locally well posed in for s > 2. The main reason is that under these assumptions 
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the operator diTij is anti-symmetric, and we are able to write 2T = C, using the notation of (3.3) and (3.4). 
The commutator term C would enable one estimates to close at the level. 

Lastly, we mention that for analytic initial data, the non-diffusive MG equation is indeed locally well- 
posed in the class of real-analytic functions, since each term in the equation loses at most one derivative. 
This was proven in [FVl lb. Theorem 2.1] by using a Cauchy-Kowalewski-type argument, in the spirit of the 
Gevrey-class regularity proof of [FT89] (see also [L097, OTOO, KV09]). In particular, this result justifies 
the short time numerical simulations of the MG equation are justified, as they only keep track of finitely 
many Fourier modes (hence the associated functions are real analytic). 

3.2. Hadamard ill-posedness for the non-diffusive MG equation. A Cauchy problem for a certain partial 
differential equation is called well-posed in X in the sense of Hadamard, if for any initial data in X, the 
problem has a unique solution in (0, T; X) , with T depending only on the initial data in the X-norm, and 
the solution map Y i— >■ L°°(0, T; X) satisfies some continuity properties, for a suitable space Y d X.\f one 
of these properties fail, the Cauchy problem is called ill-posed. In this section we discuss the failure of the 
solution map for the non-diffusive MG equation to be Lipschitz continuous with respect to perturbations in 
the initial data (in the topology of a Sobolev space X) around a specific steady profile Qq. Such results have 
been recently obtained for several models in fluid dynamics, see e.g. [GreOO, BP08, Ren09, CSIO, GVDIO, 
GVNIO, GNIO, DQSl 1, FGSV12, FR12] and references therein. 

The proof consists of a linear and a nonlinear step. We construct a steady state and a sequence of eigen- 
functions, with arbitrarily large eigenvalues, for the linearized MG equation around this steady state. Once 
these eigenvalues are exhibited, one may use a fairly robust argument (see, e.g. [Ren09]) to show that this 
severe linear ill-posedness implies the Lipschitz ill-posedness for the nonlinear problem. We emphasize that 
the key feature of the MG constitutive law (1. 12) which gives this ill-posednes is the evenness and unbound- 
edness of the symbol associated to M. In fact, the evenness also gives non-uniqueness for for bounded weak 
solutions of the non-diffusive MG equation, as was recently established by Shvydkoy [Shvl 1]. 

3.2.1. Linear ill-posedness. Let m > 1 be a positive integer. We use the steady state 

00 = sin(mx3), (3.5) 

with associated velocity field Uq = 0. Similarly to (2.22), we consider the associated Unearized (non- 
diffusive) MG operator defined by 

LO = -M3[6']93Go (3.6) 
where the M3 operator is defined in terms of its Fourier symbol (1.15)-(1.16). With this notation, we have 

Theorem 3.1 ([FVl lb]). The Cauchy problem associated to the linearized MG equation 

dtO = LO (3.7) 

is ill-posed in the sense of Hadamard over L^. That is, for any T, K > 0, there exits a real-analytic initial 
datum 6(0, •), such that the IVP associated to (3.7) does not have a solution 9 G (0, T; L?'), which obeys 

sup \\e{t,-)\\L2<K\\d{0,-)\\Y (3.8) 
te[o,T] 

where Y is any Sobolev space that embeds compactly in L?. 

The proof of the above theorem builds on the construction of unstable eigenvalues in the critically diffu- 
sive case, more specifically on Lemma 2.6. 

As in Section 2.2.2, cf. (2.23) we make the ansatz that 9 is given in terms of a Fourier sine series 

9{t,x) = e'^^ sm{kixi) sm{k2X2)'^^CpSm{mpx3), 

p>i 

with ki and k2 to be chosen precisely later. Note that since we are in the non-diffusive case, we need 
to set Ek = 0, and hence cf. (2.25) we have ap = a for any p > 1. The value of ap remain however 
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unchanged. Using Lemma 2.6, we obtain an eigenvalue = cr*(A:i, k'}), and an exponentially decaying 
sequence {cp}p>i, such that cf. (2.31) the bound 

"^^^ ^al'^i + ^2) 

holds. We set ki = and /c2 = j in (3.9), where j > 2m is a positive integer, and obtain the lower bound 

on the eigenvalue, associated to this choice of {ki,k2). Here Ci = Ci(m, A^). We also denote by 6^^^ 
the eigenfunction constructed in this way. Using an additional argument, which requires suitably bounding 
the exponentially decaying sequence {cp}p>i which defines our eigenfunction, one may in fact additionally 
normalize 9 and prove (cf. [FVllb, Lemma 2.7], [FRV12, Lemma 4.3]) that 

||0(^-)(O,-)||y = l, and ||^0-)(o, .)||^, > (3.11) 

where C2 = C2{m, N, Y) is a sufficiently large constant. For example, if y = i?^ with s > 0, we may 
explicitly compute a suitable such constant. 

We conclude the proof of Theorem 3. 1 by using a contradiction argument. If such a solution would exist, 
for a fixed T and K, we pick j = j{T, K, m, N, Y) large enough so that 

exp(rjy(2Ci)) ^ 

— — 

holds. The linear problem (3.7) can be shown to have a unique solution, and hence by construction we 
have 

, x) exp((T=i,ji), 

the assumed estimate (3.8) is contradicted by our choice of j and the lower bound (3.10). This concludes 
the proof of Theorem 3.1. 

3.2.2. Nonlinear ill-posedness. Before stating our main ill-posedness result, let us recall the definition of 
Lipschitz well-posedness. 

Definition 3.2 ([GNIO, FVllb]). Let Y C X C W^^^ be Sobolev spaces. The Cauchy problem for the 
non-diffusive MG equation is locally Lipschitz {X, Y) well-posed, if there exist continuous functions T, K : 
[0,oo)2 ^ (0, 00), so that for every pair of initial data 6^^^ (0, •), 0*^^^ (0, •) ^Y there exist unique solutions 
qW ^ 0(2) g L°°{0, T; X) of the initial value problem associated to (3.1)-(3.2), that additionally satisfy 

.) _ e(^\t, ■)\\x < K\\9W(0, •) - 0(2) (0, .)||^ (3.12) 

for every te[0,T], w/z^'r^' T = T(||0«(O, OI|y, ll^^'Ho, Oily) W K = i^(||0W(O, Oily, ||^(')(0, Oily)- 

We note that Defintion 3.2 allows the solution map to lose regularity. Since the MG equation is first order, 
the typical pairs of spaces {X, Y) that we have in mind here are X = H^, and Y = H^^^, with s > 1 + 3/4. 
The main result of this section is then: 

Theorem 3.3 ([FVl lb]). The non-diffusive MG equations (3. l)-(3.2), are locally Lipschitz {X, Y) ill-posed 
in Sobolev spaces Y C X embedded in W^'^, in the sense of Definition 3.2. 

The proof of the above result is based on the severe linear ill-posedness obtained in Theroem 3.1, and a 
standard perturbative argument [Ren09, Theorem 2] or [Tao06, pp. 183], which we breifly sketch below. 

Let ©0 = sin(mx3) be the steady state defined above in (3.5), let L6 = — M3[6']9r0o be its associated 
linearized MG operator, and denote the MG nonlinearity by N{6) = Mj[6]dj9. By our assumptions on X 
and Y , we have 

\\Le\\L2<Ci\\e\\x and \\N{e)\\L2 < Ci\\e\\\ (3.13) 



14 



SUSAN FRIEDLANDER, WALTER RUSIN, AND VLAD VICOL 



for some constant Ci = Ci(0o) > 0. Also, let Co = ||Go||y. 

In definition (3.2), we chose 6'(2)(o) = Go = e^^\t) be the steady state, and let e^'^\^) = ef\Qi) = 
Oo + SipQ, for some smooth function V'o to be chosen later. Assume by contradiction that the problem is 
Lipschitz locally well-posed in {X, Y). One may then show that there exist T = T{Co) and K = K{Co) 
such that 

sup We'iHt) - QoWx < K5 
ie[o,T] 

for any 6 G (0,5o(Co)]. Writing the solution 0^'^\t) as an order 6 perturbation around the steady state, 
6i(i)(t) = Go + 5ipsit) we thus have 

sup \\Mt)\\x<K (3.14) 
te[o,T] 

and Tps solves the Cauchy problem 

dti^s = Li^s + SNiijs), MO) = V'o- (3.15) 
The key observation is that when combined with (3.13), estimate (3.14) yields 

d\\N{i;s)\\L^<Ci6K^ 

for any t E [0, T]. Therefore, when passing 5 — )• 0, the nonlinearity in (3.15) disappears, and we should ob- 
tain the linear limiting equation. Indeed, from (3.14), (3.15), and the bounds (3.13), we obtain the necessary 
estimates to apply the Aubin-Lions compactness argument and show that tps ^ strongly in L°°{0, T\ L^) 
and on this time interval the limiting function obeys 

dt^ = Lij, V(0) = V'o. (3.16) 

But ip also inherits the bound (3.14) from the sequence ips, which contradicts Thereom 3.1 - choose ipQ an 
eigenfunction of L as constructed in Section 3.2.1, so that in time T/2 the exponential growth of the 
norm caused the solution to exceed the level 2K. This concludes our sketch of the proof of Theorem 3.3. 



4. Fractionally diffusive MG equation and well-prepared initial data 

The fractionally diffusive MG equation, although non-physical, presents a very challenging and interest- 
ing mathematical problem. Replacing the Laplacian in (1.11) by —(—A)''', for 7 G (0, 1), we obtain the 
active scalar equation 

dte + u-ve = e^{-^fe + s (4.1) 

u = M[e], (4.2) 
0(0, •) = 00 (4.3) 

where the symbol of the Fourier multiplier M is given by (1.13)-(1.16). Note that here we work with the 
non-dimensional variables, so that the initial data and the source term are normalized to be of size unity. 

The main feature of this situation is that the structure of Mj (anisotropy, evenness, and unboundedness) 
and a relatively weak smoothing effect of the nonlocal operator (—A)''' produce a sharp dichotomy across 
the value 7 = 1/2, as follows. Ifl/2<7<1 the equation is locally well-posed, while it is Hadamard 
Lipschitz ill-posed for < 7 < 1/2. At the threshold value 7 = 1/2 the problem is globally well-posed 
for 3> 1, respectively ill-posed for ^ 1- Furthermore, the anisotropy of the symbols Mj can be used 
to improve the regularity of the solutions when the initial data and source have a suitably chosen Fourier 
support. For such well-prepared data the local existence and uniqueness of solutions can be obtained for all 
values 7 G (0,1), and the global existence holds for all initial data when 7 G (1/2,1). We remark that in all 
above mentioned regimes, the problem is super-critical in the sense described in Section 2. 
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4.1. Regime of well-posedness in Sobolev spaces. In Section 2 we stated the well-posedness result when 
the diffusion in the equation is given by the full Laplacian. If the diffusion in the equation is introduced by 
the fractional operator (— A)'^, 7 G (1/2, 1), then we have the following result. 

Tiieorem 4.1 ([FRV12]). Let 7 G (1/2, 1), and fix s > 5/2 + (1 - 27). Assume that 9o G H'iT'-^) and 
S G L°^(0, 00; H^~'^ {T"^)) have zero-mean on T'^. Then there exists a time T > and a unique smooth 
solution 9 G L°°(0, T; H'{T^)) D L^{0, T; H'+^'{T^)) of the Cauchy problem (4.1)-(4.2). 

Under an additional assumption on the size of (equivalently, the size of the initial data in the unsealed 
equations), we obtain a global existence result. 

Theorem 4.2 ([FRV12]). Let 7 and S be as in the statement of Theorem 4.1, and let 9q G H^{T^) have zero- 
mean on Y^, where s > 5/2 + (1 — 27). There exists a large enough number C, such that if Sk > C, then the 
unique smooth solution 9 of the Cauchy problem (4.1)-(4.2) is global in time, i.e. 9 G L°°(0, 00; H^(T'^)). 

Remark 4.3. Note that in [FRV12] we did not use the non-dimensional variables with which we are working 
here. Therefore, the global well-posedness condition stated in [FRV12] is that the initial data and source, 
have to be sufficiently small, compared to the diffusivity coefficient k. In this paper however, the data and 
source have been normalized to be of order unity, and the only coefficient remaining is e^- It is hence 
natural that the small data condition turns to a large enough condition. This remains valid throughout 
this section. 

The detailed proofs of Theorems 4.1 and 4.2 are given in [FRV12], and rely mainly on energy methods 
and the use of suitable commutator estimates. We omit these details here. 

4.2. Regime of ill-posedness in Sobolev spaces. In this subsection we consider the stability of steady 
states to (4. l)-(4.2). The construction follows the outline described in Section 2.2.2. In this case, we choose 
00 = sin(mx3) and S = EkIti?'^ sm{mxs), for some integer m > 1. As before, the first step in proving 
that the Cauchy problem (4.1)-(4.2) is Lipschitz ill-posed is to consider the linearized problem around the 
given steady state, and exhibit the existence of unstable eigenvalues with arbitrarily large real part. Thus, 
we consider 

dt9 = L9, (4.4) 

where the Unear operator L is given by 

L9{x,t) = -M:i9{x,t)d39o{x3) - e^{-Ay9. 
We seek the temperature 9 which solves (4.4) and has the form 

9{t, x) = e"^ sin(fcixi) sin(A:2X2) Cp sm.{mpx3) (4.5) 

p>i 

for some fixed ki,k2 G to be chosen later. Inserting (4.5) into (4.4) and using the definition of M^, we 
obtain that (2.24) holds, but with cxp denoted by 

ap = a + e^{kl + kl + {mpfy (4.6) 

and with ap given by (2.26) above. As before, these assumptions lead to solving the continued fraction 
equation (2.29) for some real o" > 0. This is achieved by suitably modifying Lemma 2.6 to obtain: 

Lemma 4.4. Let ap, ap be defined by (2.26) and (4.6) respectively, and let ci = ai. There exists a real 
eigenvalue cr* > 0, and a rapidly decaying sequence {cp}p>2 so that the function 9 defined by (4.5) solves 
(4.4). Furthermore, we have the lower bound 
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In particular, setting for instance ki = and k2 = j in (4.7), for each j > 2m we obtain a real eigenvalue 
f7=K J of L, which is bounded from below as 

Moreover, given a Sobolev space Y that is compactly embedded in L^, we may ensure the existence of a 
constant C2 = C2{m, N, Y, 7) such that the smooth eigenfunction 9^^^ corresponding to cj* j obeys (3.11). 
The immediate consequence of estimate (4.8) is that due to 1 — 27 > 0, we have 

(T=Kj — )• 00, as j — )• 00. 

The above construction proves that the linearized problem (4.4) is ill-posed in L^, as we have exhibited 
eigenvalues of arbitrarily large real part. Following the argument outlined in Section 3.2.2, this type of 
strong ill-posedness for the linearized problem, automatically implies that the nonlinear problem is Lipschitz 
ill-posed in Sobolev spaces (recall Defintion 3.2 of Lipschitz well-posedness). Thus we obtain: 

Theorem 4.5 ([FRV12]). The Cauchy problem (4.1)-(4.3), with 7 G (0, 1/2) is locally Lipschitz. {X,Y) 
ill-posed in Sobolev spaces Y <Z X embedded in M^^'^(T'^). 

4.3. The threshold case 7 = 1/2. The value 7 = 1/2 marks the threshold for proving the local well- 
posedness in Sobolev spaces for (4.1)-(4.3). If the initial datum is small with respect to k, or equivalently, if 
is sufficiently large (see Remark 4.3 above), then one may use energy estimates to show that there exists 
a unique global smooth solution evolving from this data. In this direction we have: 

Theorem 4.6 ([FRV12]). Let s > 5/2 and assume that 6q G H'^iT^) and S G L;^(0, 00; H''-^{T^)) have 
zero-mean on T^. There exists a sufficiently large number C, such that if > C, then the unique smooth 
solution 9 G L°°(0, T; H^(T'^)) of the Cauchy problem (4.1)-(4.3) is global in time. 

As in the case of 7 G (1/2, 1), the proof is based on the classical Kato-Ponce commutator estimate. The 
proof however does not apply for the case of "small" (large data), not even to prove the local existence 
of solutions. 

In the case where is suitably small, we may construct a steady state such that solutions evolving from 
initial data which is arbitrarily close to this steady state diverge from it at an arbitrarily large rate, i.e. the 
problem is Hadamard ill-posed. More precisely, we have: 



Theorem 4.7 ([FRV12]). Let Y C X C ^^^'''(T^) be Sobolev spaces, let 7 = 1/2, and fix the integer 

m = 1. Consider the steady state s< 
S = eK?Tisin(mx3), and assume that 



171 = 1. Consider the steady state solution of (4.1) given by 9^^ = sin(mx3) G Y, with source term 



(2) 

Then, given any T > and any K > 0, there exists an initial condition 9q ^ Y and a corresponding 
solution 6'^^) G L°°(0, T; X) of the Cauchy problem (4.1)-(4.2), such that 

sup \\9(^\t, •) - 9W(t, OIU > 2K||ef - 9^^^\\y 
te[o,T] 

where 9^^\t, •) = ^g^^. That is, the problem (4.1)-(4.2) is locally Lipschitz (X,Y) ill-posed when ^ 1- 

To prove the above result, we follow the ill-posedness result in the case 7 G (0,1/2). The main difference 
is that while for 7 G (0, 1/2) by (4.8) we have real eigenvalues of arbitrarily large real part, when 7 = 1/2 
the lower bound (4.8) only blows up for sufficiently small, depending on N'^ and our choice of m. This 
condition is rather unsatisfactory, since it implies — > as A^^ — )• 0. In fact, we can estabUsh a more robust 
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smallness condition on in the sense that it is independent of N'^ (and m). To see this, insert ki = j'^N ^ 
and k2 = jm^^"^ in the lower bound (4.7), where j > max{A^^m, 4m}. We obtain the lower bound 

Therefore, under the assumption (4.9), which is independent of m or N, the right side of (4.10) diverges 
to 00 as j — )• 00 and hence we again have a sequence of unstable eigenvalues of the linear problem, with 
arbitrarily large real part. Combined with the nonlinear argument presented in Section 3.2.2 this concludes 
the sketch of the proof of Theorem 4.7. 

4.4. Well-prepared initial data. In this subsection we present results which take advantage of the curved 
nature (paraboloid-type) of the regions in frequency space where the symbol M{k) grows like We 
observe that if the support of 9 lies on a plane, then the operator M behaves like an operator of order zero 
acting on this 6, and hence the velocity u is as smooth as 9. This allows us to improve the well-posedness 
results in the range > and < 7 < 1, over the generic case when no conditions on the Fourier 
transforms of the initial data (and source term) are imposed. 

To be more precise let us introduce some notation. For a rational number g G Q with g / 0, define 
frequency plane Vq as 

Vq = {k = (fci, fe2, A;3) el? ■.k2 = qki}. 

One may show that whenever g / 0, 

\Mj{k)\<Cq< 00, for all /c G and j G {1,2,3}. (4.11) 

This observation, combined with the fact that when the data and source S have frequency support on Vq, 
then so does the solution 9{t, ■) of the Cauchy problem (4. l)-(4.3) at all later times i > 0, means that the 
smoothing properties of the fractional diffusion term are stronger than they are for the generic data situation. 
Hence it is possible to prove stronger regularity results, and in particular the local existence and uniqueness 
of Sobolev solutions holds /or all values of 7 G (0, 1). 

Theorem 4.8 ([FRV12]). Let 7 G (0,1), and fix s > 5/2 - 7. Assume that 9o G H%T^) and S G 
L°°(0, 00; H^~'^ (T"^)) have zero-mean on T'^ and are such that supp(0o) C Vq and supp(S'(t, •)) C Vq 
for some fixed q G Q \ {0} and all t > 0. Then there exists a T > and a unique smooth solution 9 G 
L~(0, T; F*(T3)) n L2(0, T; H-'+^{T^)) of the Cauchy problem (4.1)-(4.3), such that supp(^(t, •)) C Vq 
for all t£ [0,T). 

We note that the above theorem holds even if7G (0,1/2), a setting in which we know that for generic 
initial data the problem is ill-posed. Moreover, one may also prove that if the initial data is small with respect 
to K, or equivalently, » 1 in the non-dimensional setting of this manuscript, the solution obtained in 
Theorem 4.8 is global in time. 

The main difficulty in the proof of Theorem 4.8 is the construction and uniqueness of solutions, which live 
on the frequency plane Vq. This relies on constructing of a suitable sequence of approximations supported 
on the Fourier space on Vq and showing its convergence using energy methods. The main obstacle one has 
to overcome is designing an iteration scheme which is both suitable for energy estimates, and is such that 
in each iteration step the frequency support of the approximation lies on Vq. More explicitly, a scheme 
in which dt9n+i is given in terms of v ■ V9n preserves the frequency support on Vq, however is not good 
for closing the estimates at the level of Sobolev spaces. On the other hand if we define dt9n+i in terms of 
v-'V9n+i, the energy estimates are now clear but it seems difficult to inductively prove that 9n+i is supported 
on Vq. Therefore, we work with a regularized problem, where we introduce an artificial term —eA9. The 
advantage of considering a stronger dissipative problem is that we can construct smooth solutions which 
have frequency support lying on Vq. For details of the construction we refer the interested reader to [FRV12]. 
Lastly, we pass to the Umit e — ;> and obtain a solution of (4.1)-(4.3) with desired properties. 
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If the initial data 9q and the source term S have Fourier support in Vq, we show that for 7 G [1/2, 1) the 
local in time solution may be continued for all time. 

Theorem 4.9 ([FRV12]). Let s,Oq and S be as in the statement of Theorem 4.8. If ^ > 1/2, the unique 
smooth solution 6 of the Cauchy problem (4.1)-(4.3) is global in time. 

The proof hinges on the observation that the solution constructed in Theorem 4.8 satisfies supp(^) C Vq. 
Then, in view of (4.11) and the Hormander-Mikhlin theorem we have HuHh/s p < C||^||vi/s,p for all 2 < 
p < 00, s > 0, and so the velocity u is as smooth as the scalar 9. Hence, the case 7 > 1/2 is "sub-critical" 
for the a priori controlled quantities, and energy estimates in the spirit of those used to prove the global 
regularity of the sub-critical SQG equation [CW99, CC04, Wu05], may be used to prove Theorem 4.9. The 
case 7 = 1/2 in Theorem 4.9 is more subtle, and is analogous with the problem of global regularity for 
the critical SQG equation [CVIO, KNV07, KN09, CV12]. The anisotropic nature of the symbol M, makes 
it more convenient to apply the De Giorgi-inspired proof of [CVIO] in the context of the fractional MG 
equation with 7 = 1/2. The main ingredient is that if 6 € L°° and supp(0) G Vq, then u G BMO. Hence, 
considering first the linear problem with divergence-free BMO drift we obtain Holder continuity of the 
solution. The later norm is sub-critical and one may bootstrap to higher regularity. We refer the interested 
reader to [FRV12] for details. 

5. Conclusions 

We have proven that the magneto-geostrophic model proposed by Moffatt and Loper is capable of sus- 
taining dynamo action. The MG-dynamo instability exhibited here relies crucially on the fact that the di- 
mensionless parameter (1.5) is neither 0, nor 00. The exponential growth is of order 0{exp{t/eK)), for 
all A^^ G (0, 00), and < ^ 1. In terms of the physics in the model, this condition on N"^ requires that 
all three forces, Coriolis, Lorentz, and gravity, are non-zero, making the MG-instability somewhat different 
in nature from the usual Rayleigh-Taylor instability. 

At the level of mathematical analysis several problems remain to be addressed, including some well/ill- 
posedness issues. For the non-diffusive equation, is there a "stronger" form of ill-posedness such as norm 
inflation or even finite time blowup? Can we refine the well-posedness results for the fractionally diffusive 
problem in the range 1/2 < 7 < 1? That is, can one establish the global existence of smooth solutions or 
on the contrary, prove finite time blowup? 

Several, interesting, mathematical issues arise when returning to the derivation of the MG equation from 
the full three dimensional MHD equations in the spirit of [ML94, Mof08]. In particular, we note that 
Moffatt's and Loper's initial discussion includes the kinematic viscosity v, which is directly proportional 
to £ in (1.1). The values of e and are highly conjectural, but both are believed to be extremely small. 
Therefore, as opposed to the MG system, if we keep e > we derive from (1.1)-(1.4) the system 

Ar2e3 X n = - VP + 62 • V6 + 7V^(9e3 + eAu (5.1) 

= 62 • Vn + A6 (5.2) 

dte + u-Ve = e^M + S (5.3) 

V-n = 0,V-& = 0. (5.4) 

The presence of the Laplacian term in (5.1) makes the constitutive law relating u and 9 "better", in the sense 
that the M operator corresponding to (5.1)-(5.4), now depends both on N and e, and is in fact smoothing 
of degree 2 (see equations (29)-(30) in [ML94]). Hence, in this case the constitutive law is even better than 
that for the two-dimensional Euler equation in vorticity form, implying that: in the range e > 0, the system 
(5.1)-(5.4) is globally well-posed, even if £k, = 0. 

This observation is interesting in view of the much more challenging problem of establishing properties 
of the solutions to (5. l)-(5.4) in the vanishing parameters limits e — )■ and 0, as proposed by Moffatt 
and Loper. The later question is directly related to the turbulent scaling regime for 9 conjectured by Moffatt 
in [MofOS]. "Turbulence" in the framework of vanishing diffusivity for the forced system (5.1)-(5.4), 
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and the presence of anomalous dissipation of energy for 9, or "thetargy", may be plausible for very weak 
solutions (below Onsager critical) of the non-diffusive equation. 
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